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Q ; Abstract 
cn 

We study solutions of statistically neutral polyampholyte chains containing a large fraction of neu- 
CN . tral monomers. It is known that, even if the quality of the solvent with respect to the neutral 

, monomers is good, a long chain will collapse into a globule. For weakly charged chains, the interior 

Q^ \ of this globule is semi-dilute. This paper considers mainly 6*— solvents, and we calculate the dy- 

^ I namic charge density correlation function g{k, t) in the interior of the globules, using the quadratic 

approximation to the Martin-Siggia-Rose generating functional. It is convenient to express the 
results in terms of dimensionless space and time variables. Let ^ be the blob size, and let r be the 
characteristic time scale at the blob level. Define the dimensionless wave vector q = S,k, and the 
\ dimensionless time s = tjr. We find that for q <\, corresponding to length scales larger than the 

blob size, the charge density fluctuations relax according to g{q,s) ~ q^{l — s^/^) at short times 
• s < 1, and according to g{q,s) ~ at intermediate times 1 < s < We expect these 

^ ■ results to be valid for wave vectors q > 0.1, where entanglements are unimportant. 

O 
> 

X ■ I. INTRODUCTION 

A polyampholyte is a polymer chain that contains electrically charged monomers of both 
signs. Apart from the charged monomers, the chain may also contain neutral monomers. In 
this paper we study polyampholytes in which the neutral monomers are in the majority, so 
the chains are only weakly charged. The positive and the negative monomers are distributed 
with equal probability and without correlations along the chain, separated from each other 
by flexible neutral spacers. The chains are present in a solvent. We consider two situations; 
either the quality of the solvent with respect to the neutral monomers is good, or the solvent 
is under ^—conditions. It was shown in |1|] that in both situations an isolated neutral chain 
collapses into a non-compact globule, provided that the chain length exceeds some critical 
value. For weakly charged chains the polymer concentration inside the globule is low. In 
case of a good solvent, further collapse is prevented by the second order virial coefficient 
(also called the excluded volume parameter), whereas in case of a 6*— solvent, further collapse 
is prevented by the third order virial coefficient The non-compact semi-dilute globule 
can be regarded as a dense melt of blobs. At length scales small compared to the blob 
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size ^, the electrostatic interaction is weak (relative to the entropy), and the system is 
indistinguishable from a neutral non-collapsed polymer coil. In case of a good solvent, the 
chain conformation is a self-avoiding walk, whereas in case of a ^—solvent, it is a random 
walk. On these short length scales, the charged monomers occupy random positions in 
space, and so the electrostatic interaction is unscreened. The importance of the electrostatic 
interaction relative to the entropy increases if one considers larger and larger length scales, 
and the break-even point is at the blob size. At length scales larger than the blob size 
the interaction is strong, and the charged monomers rearrange themselves spatially in order 
to minimize the energy. This leads to a screening of the interaction. A side effect of these 
rearrangements of the charges is an effective attraction, causing the chain to appear collapsed 
on length scales larger than ^. This physical picture, originating from Ref. [H], shows why 
the electrostatic screening length has to be equal to the blob size ^, which is also equal 
to the screening length of the excluded volume interactions in case of a good solvent. Since 
for weakly charged polyampholytes the blob size is much larger than the average distance 
between two charged monomers, Debye-Hiickel theory is applicable (see also Ref. If 
the chain has a net charge, there is the possibility that the formation of the globular state is 
prevented by the electrostatic repulsion, depending on the chain length and the magnitude 
of the excess charge. For the "statistically neutral" chains considered here, the typical excess 
charge is proportional to the square root of the chain length. It was shown in Ref. that 
in this case the electrostatic repulsion is insufficient to prevent the formation of the globular 
state for long enough chains. The only effect of the presence of a relatively large excess 
charge is that the globule becomes elongated. In Ref. [|] multi-chain effects were studied. It 
was shown that solutions with non-zero concentration start to phase separate already at very 
low concentrations, even if the quality of the solvent with respect to the neutral monomers is 
good. The supernatant consists of isolated spherical globules floating in the solvent. These 
globules contain either a single neutral chain, or two chains whose excess charges cancel each 
other. The semi-dilute precipitate, which still has a low polymer concentration, is locally 
indistinguishable from the interior of the globules. In this paper we will concentrate on the 
dense (semi-dilute) phase. 



II. METHOD AND REGION OF APPLICABILITY 

We study the dynamics of the charge density fluctuations in the dense phase (precipitate 
or globule), using the quadratic approximation to the Lagrange version [5,6] of the Martin- 
Siggia-Rose formalism [0 . We will follow closely the techniques of a paper by Fredrickson and 
Helfand [Q. Readers who wish to follows the details of our calculation are advised to study 
Ref. 1^ first, especially its appendix. First we will determine in what length scale regime 
this formalism is expected to give correct results. There are three important length scales 
in the system: the blob size ^, the hydrodynamic screening length ^h, and the entanglement 
length ^e- At length scales smaller than ^, the thermal energy is larger than the electrostatic 
energy, and the system behaves as if it were a neutral single chain. Since in the past the 
dynamics of neutral chains has been studied extensively 0, and we are mainly interested 
in the influence of the electrostatic interaction, we will only consider length scales that are 
larger than C,- When observed at these scales, the chain conformation is always a random 
walk, even if the quality of the solvent is good. This is due to a screening of the excluded 
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volume interactions 0. It means that the excluded volume parameter v enters the physics 
at large length scales > ^ only by fixing the blob size 

Another important length scale is the entanglement length ^g- It is very difficult to 
take the inffuence of entanglements on the dynamics into account, and we will not make 
an attempt to do so. This means that our results will not be valid on length scales large 
compared to the entanglement length ^g- One can obtain an estimate for in the following 
way. In a polymer melt, the number of monomers A^^e in between two entanglements is in 
the range A^^g ~ 50 — 500. Since the interior of the polyampholyte globules can be regarded 
as a dense melt of blobs, it follows that the typical entanglement length is given by 

^g ~ ^Ni'^ ~ loe (1) 

At length scales that are large compared to it is to be expected that the dynamics 
occurs via some kind of reptation, and on these scales the relaxation of charge density 



ffuctuations probably becomes very slow (compare with Ref . |T0[ , where systems with short 
range interactions are studied). 

Next we discuss the hydrodynamic screening length ^h- At length scales smaller than ^h, 
the chain exhibits Zimm-like behavior, whereas at length scales larger than ^h, it exhibits 
Rouse-like behavior. Since it was shown in Ref. 0] that the quadratic approximation leads to 
the prediction of Rouse-like behavior, our analysis is only applicable to length scales larger 
than ^H- We need an estimate of its value. According to Ref. 

«H~^, (2) 

where p is the monomer number density, and h is the statistical segment length of a 
monomer (NB in this paper the sign ~ means: has the same order of magnitude). First 
consider the case of a ^—solvent. According to Ref. |l|] we have 

where £ := / ATiek-oT is the Bjerrum length (£ = 0.7nm in water at room temperature), 
e is the electric permittivity of water, k-B is Boltzmann's constant, T is the temperature, and 
/ is the fraction of monomers that carries a charge. Combining Eqs. (1) and (2), we obtain 

^H,6l ~ 6- (4) 
In a good solvent, characterized by excluded volume parameter v, we have Jll 

Combining Eqs (1), (2) and (4), we obtain 



en, 
4 




(6) 



3 



We conclude that under good solvent conditions the hydrodynamic screening length is 
larger than the blob size, whereas under 6'— conditions they are approximately equal. 

Finally, we discuss the region of validity of the quadratic approximation. It entails 
expanding the effective Hamiltonian in powers of the charge density ip{r), and retaining only 
the second order term. This approximation can only be justified if the typical amplitude of 
ip{r) is small. At length scales smaller than ^, the system is strongly fluctuating, and the 
amplitude of iplr) is large. At length scales larger than ^, the charge density fluctuations 
are strongly suppressed by the electrostatic interactions, and the amplitude of ip{r) is small. 
It follows that the quadratic approximation is acceptable at length scales larger than ^, but 
not at length scales shorter than ^. 

Summarizing, the approximations made in this paper are expected to be reasonable on 
length scales that are large compared to both ^ and ^h, and small compared to ^e- In order 
to ensure that our region of applicability is not empty, we will henceforth assume that the 
solvent is under ^—conditions, in which case ~ ^, and ~ 10^. 

Since in a 6'— solution the hydrodynamic interaction is screened beyond length scale ^, 
its influence enters the physics on intermediate length scales ^ < L < in a trivial way via 
just one parameter, which is the characteristic time scale r at the blob level. Since inside 
the blobs the influence of the electrostatic interactions is negligible and the chain is not 
collapsed, r can be found by considering a single neutral ideal chain with hydrodynamic 
interaction (Zimm- model). This leads to the estimate 

where t] is the viscosity of water. 



III. CALCULATION OF THE DYNAMIC CHARGE CORRELATION FUNCTION 

As mentioned in the introduction, we wish to describe the (semi-dilute) precipitate of a 
solution of weakly charged polyampholytes. The individual chains are assumed to be much 
longer than the minimum length required for an isolated chain to collapse into a globule. 
All chains have exactly the same number of charged monomers, and these monomers are 
placed at regular distances along the backbone. The last two assumptions, which are made 
in order to simplify the model, do not restrict the generality of our results. Let N be the 
number of charged monomers per chain, e the charge per monomer. Hp the number of chains 
in the precipitate, and V the volume of the precipitate. Note that V cannot be chosen 
independently from Up and A^, because the precipitate has a well-defined density Let 
the coarse grained conformation of chain a be described by the 3-dimensional vector Ra(T) 
, where r is a continuous parameter running along the backbone. It is defined such that 
for two charged monomers neighboring along the chain we have At = 1. We simplify the 
model by smearing out the electric charges along the chain, in such a way that e6{T)dT is the 
amount of charge in between the points labeled by r and r + dr, where 6{t) is a Gaussian 
random variable with first two moments 

(eair)) = 0, {9a{T)9,{r')) = 6{t - r')5,,. (8) 
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Since we will work in the quadratic approximation, the fact that we switch to a Gaussian 
charge distribution has no effect on the final result. We will describe the system by means 
of the following effective Hamiltonian: 



n 




(9) 



where h is the root-mean-square distance between two charged monomers neighboring 
along the chain. From now on, we will choose the units of length, mass and time such that 

7 = 1, A;bT = 1, 62^3, (10) 

where 7 is the effective friction coefficient per charged monomer. Ultimately, we are 
interested in how fluctuations in the charge density decay in space and over time. In terms 
of the annealed variables Ra(T) and the quenched variables 6 air), the charge density is given 
by e^/'(r,t), whereQ 

^Pir,t)=Y: fdr9^ir)6ir-R,ir)). (11) 

In terms of ip, the Fourier transform of the dynamic charge density correlation function 
is given by 

^,(k,t) = ^(^(-k,t)^(k,t)), (12) 

where the brackets denote an average over the annealed variables. Note that the cor- 
relation function depends explicitly on the disorder. The calculation is simplified consid- 
erably by the assumption that the correlation function is self-averaging, which means that 
gg(k,t) = g{k,t) := gg'(k,t) with probability 1, where the bar denotes an average over the 
quenched variables 6'. This means that in order to find geO<; t) for a representative 6 drawn 
from the Gaussian probability distribution Eq. (^ (which is what we are after), it suffices 
to calculate the average of this quantity over the disorder. 

We will calculate 5'(k, t) by means of the Lagrange version of the MSR (Martin-Siggia- 
Rose) formalism. In the appendix we have provided a brief derivation of the MSR functional. 
It is convenient to introduce an external field h{k,t) that couples to '?/'(k, t). The dynamic 
charge correlation function can then be found by differentiating the logarithm of the MSR 
functional Z[h] twice with respect to h. In principle, the average over the quenched disorder 
should be performed after the differentiation, but since Z[h = 0] = 1 is independent of the 
quenched variables (see the appendix), the order of operations can be interchanged. This is 



^The hat on i/j denotes the dependence on the annealed variables Ra(r), and the quenched vari- 
ables 9a{T). This notation should not be confused with the notation in the MSR- formalism, where 
the hat indicates conjugated variables (see appendix A). We stick to these notations because of 
convention. 
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a great simplification, and it makes the MSR formalism especially suitable for systems with 
quenched disorder. 

Starting from Eqs. ( |A9| ) and ( |A10| ), and substituting the effective Hamiltonian Eq. (|[) 
for the interaction energy , one arrives at the expression for the MSR functional. Since 
analogous calculations have been published before (see in particular Ref. we will not 
present the details, but just give the result. During the derivation, the following fields arise: 

^i(k,t) = ^ fdT9a{r)exp [2k-R,(r,t)] 

V^2(k,t) = 5] /dr^^,(r)k-R„(r,t)exp[zk-R,(r,t)]. (13) 
In the quadratic approximation, the final disorder averaged expression for Z[h] is 

Z[h]oQ J VipiVij2'D(j)iV(j)2ew + J J^h{-k, -u) ^i{k, u) , (14) 
where the Lagrangian C is given by (there are summations over the indices i,j) 
C = -ij j ^i)i{-k,-uj)^2%^) -i f f (l)i{-k,-uj)ipiik,uj) 

J u) rC Jul Jk 

+^ / / V^j(-k,-a;)^i(-k,-^)^,(k,^). (15) 

Z Jllj Jk 

c = UpN/V is the density of charged monomers. The integral measures are defined by 
= i/d- . (16) 



2tt J 7k ■ (27r)3 . 

Although in the full expression for the functional Z[h] there are more fields present, in 
the quadratic approximation these fields couple neither to ipi, nor to (pi, and can there- 
fore be omitted. One consequence of this is that in this approximation the hydrodynamic 
interactions do not have an influence on the dynamic charge density correlation function. 
Therefore, our results are only valid on length scales larger than the hydrodynamic screen- 
ing length. The integrals over ipi and (pi in Eq. ( JT^ j are Gaussian and can be calculated 
explicitly, after which g{k, t) can be obtained by differentiating the result twice with respect 
to the external field h. The result is 

/, N cVu{k,uj) 

'""^ {l-^i''k-^Vu{k,u;)){l-K-'k-^V2i{k,u)y ^ ^ 

where := Airic is the Debye-Hiickel expression for the inverse square screening 
length. The functions Vij are given by 

Vii{k, t) = 1 1 dr (exp [zk ■ (R(r, t) - R(r, O))])^ 

Vuik, t) = 1 1 dr (k ■ R(r, t)exp [zk ■ (R(r, t) - R(r, 0))])^ 

^21 (k, t) = -^l dr (k • R(r, 0)exp [zk ■ (R(r, t) - R(r, 0))])^ 

V22{k, t) = -j^J dr (k • R(r, t)k ■ R(r, 0)exp [^k ■ (R(r, t) - R(r, 0))])^ . 

(18) 
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The average (■ 
which is given by 



■)g is calculated with respect to the single-chain Rouse Lagrangian Cn 



C = I I drR(r, -uo) ■ R(r, u) + f f d™R(r, 



-UJ 



R(r,^) 



drR(r, 



-UJ 



(19) 



Before starting to calculate Vij(k, t), it is useful to determine what are the relevant length 
and time scales. As discussed in the introduction, we are interested in the processes occurring 
at length scales of the order of the blob size ^. As is usual for the semi-dilute regime, the 
blob size reaches a finite limit when the chain length approaches infinity. This limiting value 
is reached once the chain length exceeds the value necessary for an isolated chain to collapse 
into a globule. It follows that for the processes occurring at length scale ^ the limit N ^ oo 
is meaningful, and approached easily in experimental situations. Therefore, we can safely 
assume that the wave vectors of interest satisfy k ~ >> N~2, which will simplify the 
calculations. Considering the time scale, it is clear that the processes occurring at length 
scale ^ are much faster than those occurring at length scale N2 . Therefore, in the calculation 
of Vi,j{k, t) we can restrict ourselves to times satisfying t « t^, where is the Rouse time 
0] of a single chain. In this wave vector and time regime it is possible to find explicit and 
simple expressions for Vij(k, t). Those readers who are interested in the calculation leading 
to Eq. (|20|) are advised to study the appendix of Ref. |@], where similar calculations are 
worked out in detail. The result is 



Vu{k,t) = exp 
Vu{k,t) = -e{ 

V2i{k,t) = -e{t)- 
^22(fc,t) =0, 



k' 



t|l/2 



exp 



-A;2|t|V2 
fc2|t|V2 



(20) 



where 6{t) is the Heaviside step function. We omitted numerical constants of order unity. 
Although it is possible to find explicit expressions for the Fourier transforms Vij{k,uj), the 
resulting expression for g{k, uj) would be rather complicated. Instead, it is much more useful 
to derive a transparent, albeit approximate, expression for g{k,uj). An additional advantage 
is that this simplified expression can easily be inverse Fourier transformed with respect to 
the frequency. Consider the frequency regime uj » k'^. Further on we will see that the 
characteristic frequency at wave vector k satisfies this criterion. By substituting z = —iujt, 
distorting the integration contour back to the real axis, and expanding the integrand in 
powers of zk^/uj P], one arrives at the following approximate expressions for Vij{k,uj): 



\uj\ » k 



Vu{k,uj 
Vi2ik,uj 




uj>0 
UJ <0. 



(21) 



Combining Eqs. (|T^ and ( l2l|) one obtains 
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q(k,uj)(X r ^ T \uj\»k'^. (22) 

£|CU|V2 [(1+ |cu/t-4|l/2)2^l] 

There are two frequency regimes: 

As discussed before, the formahsm used can only be expected to describe the system 



correctly at length scales larger than the blob size. Therefore, we will consider Eq. (p3D 
only for k < k. In that case, both frequency regimes given in Eq. (|2^) lie within the range 
of validity Ic^l >> k'^. Doing an inverse Fourier transform with respect to the frequency we 
find two regimes: 

c 'g{k,t) ^k\ ,4,-1/2 ,.-4 / , / ,.-4 (24) 



<t<k-^ 

with the restriction C, < k~^ < C,c on the wave vector. Note that, in accordance with 
the results of Ref. the wave vector dependence of the static correlation function g[k, t = 
0) coincides with the Debye-Hiickel expression g oc /c^/ {k"^ + k^) in the regime k << k. 
Equation ( pH ) predicts that charge density fluctuations decay with a wave vector independent 
initial rate r ~ ~ The exponent 4 is characteristic for the Rouse model. Since at 
length scales smaller than ^ the hydrodynamic interactions are not screened, one would 
expect a rate r ~ S,^. The discrepancy is due to the Gaussian approximation This does 
not mean that Eq. ( p4D breaks down completely: since at length scales larger than ^ the 
hydrodynamic interactions are screened and the chain is Rouse-like, Eq. (|^) is valid for 
k < K, provided that one inserts the correct basic time scale r at the blob level; see the 
discussion at the end of the introduction. Defining the dimensionless wave vector q := k/n 
and the dimensionless time s := t/r, we obtain for 0.1 < g < 1 

-1 / N 2(1- S <1 

c g{qK, ST)^q x | 1< s < q-^ ^ ^ 



with 



knT 



(26) 



The expression for r reflects the Zimm behavior at small length scales, whereas the 
exponent 4 appearing in the condition s < q~^ reflects the Rouse behavior at longer length 
scales. The time dependence of the dynamic charge correlation function is schematically 
depicted in Fig. 1. 



IV. DISCUSSION 

In order to determine the role of the electrostatic interactions in the expression for the 
dynamic charge density correlation function, consider a system in which the electrostatic 
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interactions between the charges are switched off. Physically, this could be achieved by 
immersing the polymer in a concentrated salt solution. Although in this case the chains 
would not collapse or phase separate, we will still assume that the system is semi-dilute, 
for instance by imposing a non-zero concentration. In this case, the charge density corre- 
lation function would follow immediately from Eqs. ( p!7D and ( pOD by taking k = . It 
follows that when the electrostatic interactions are absent, the charge density correlations 
decay according to a stretched exponential g{k,t) ~ fc^/t~^exp [-(t/r')'/'l for some T . 

1 /2 

Considering that the Taylor expansion of this function starts with 1 — {t/r') , one sees 
from Eq. p5| ) that at short times an interacting system relaxes in qualitatively the same 
way as a non- interacting system (though the decay rate is higher), but that at longer times 
the relaxation is completely different. We conclude that the appearance of the power law 
g oc must be due to the Coulomb interactions. In order to determine the influence of 

the polymeric bonds, we calculated for comparison the dynamic charge density correlation 
function of a salt solution, using the same formalism and approximations. The result is an 
exponential decay of the correlations over time: 

^) = TT^^^'P (27) 

The time dependence of Eq. ( p7| ) is completely different from that of Eq. (^Sf ), meaning 
that the presence of the polymeric bonds has a large influence on the relaxation of charge 
density fluctuations. 

The dynamic correlation function for the total density is qualitatively different from that 
for the charge density. As an illustration, consider the correlation function for a semi-dilute 
homopolymer solution, which has been calculated within the quadratic approximation in 
Ref. P]. On length scales that are small compared to the radius of gyration of the chains, 
but large compared to the correlation length ^, the fluctuations in the total density decay 
exponentially with time (see Eq. (3.19) in Ref. ||^). Experimentally, it would be interesting 
to test the existence of the two frequency regimes in Eq. (pHf) by means of scattering 
experiments. 



APPENDIX A: 

In this appendix we present a brief derivation of the Martin-Siggia-Rose functional [0. 
We will follow the method developed in Refs. PJ^. Consider a system of interacting 
point particles in solution. Let n = 1, - ■ ■ ,N number the particles, let a = 1, 2, 3 denote 
a Cartesian coordinate, and let Rna be the a— component of the position of particle n. 
Let R denote the 3A^— dimensional vector with components Rna, and let U [R] denote the 
interaction energy. In case that R represents the coarse grained conformation of a polymer 
chain, entropic contributions have to be taken into account and one should replace U by an 
effective Hamiltonian 7i. The time evolution of the probability density P [R, t] is governed 
by the Fokker-Planck equation ||9,11 



dP 

— = V ■ L ■ [fceTVP + PVU] , (Al) 
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where V represents the vector with components djdRna-, and dots denote inner products. 
The mobihty matrix L [R] is given by [0 



■'na^mfi 



HaP (Rri 



Rri 



n ^ m 



5a0 

7 



(A2) 



Hapiv) is the Oseen tensor (see Ref. 0, appendix 3. Ill) describing the hydrodynamic 
interaction, R„ denotes the position of particle n, and 7 is the friction coefficient per particle. 
In order to derive the Martin-Siggia-Rose functional it is convenient to switch first from 
the Fokker-Planck equation to the equivalent Langevin equation [|ll[]. In the Stratonovich 
interpretation it is given by p[ 



(v-g) + g-C- 



fro 1 

— = -L-Vf/ + fcBTV-L--j 

at 2 

The stochastic force is a white Gaussian noise with first two moments 

(C(t)) = o {C{t)C{t')) = i6{t-t'). 



(A3) 



(A4) 



I is the 3A^ x 3A^ identity matrix. The matrix g is related to the mobility matrix hfL 
by i 



g ■ gT = 2fcBTL. 



(A5) 



Equation (A5) leaves some freedom in the choice of g, which can be used to impose the 
condition p 

V ■ g = 0. (A6) 
Using this, the Langevin equation simplifies to (it can be shown that V ■ L = 0) 



— = -L.Vt/ + g-C + h. 

at 



(AT) 



The external force hi{t) working on particle i is introduced in order to be able to calculate 
response functions. This force will be set to zero afterwards. In Refs. [5,6] it is worked out 
in detail how one can derive, starting from a Langevin equation, the expression for the path 
probability distribution V [R(t)]. By introducing the so-called conjugate field R(t), which 
is rather straightforward in the Lagrange formalism [0], one finds the expression 



V 



R(t),h = / I?R(t)exp 



/:[R,R] + / dth(t) ■iR(t) 



where the "Lagrangian" £[R, R] is given by 



£[R,R1 



dt 



z^bTR ■ L • R + iR ■ (R + L ■ Vt/ 



(A8) 



(A9) 



The Martin-Siggia-Rose functional is defined by integrating the path probability over all 
possible evolutions R(t) of the system, in the presence of the external fields h and h, where 
the field h couples to R. One obtains 
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Z[h, h] = Jl^^J VR{t) exp -C[R, R] + y dth{t) ■ R{t) + J dth{t) ■ zR(t) . (AlO) 

It follows from Eq. ( [A8| ) that the correlation and response functions of the system de- 
scribed by the Langevin equation can be obtained from the MSR-functional Eq. (|A10 ) 



by differentiation with respect to the external fields. However, the continuum expression 



for Z, as it is given here, is ill-defined fT^. For instance, it is not possible to determine 
the value of the equal-time response function (R{t)iR{t)'^ . Retracking the derivation of Eq. 
( |A1CI|) one finds that the discretization underlying the integrals over time is such that the 
equal-time response functions are identically zero |jl2|. This extra information is sufficient 
to remove all ambuigity from Eq. (AlO). 



We discuss briefly the adequacy of the Langevin equation Eq. ( |A7| ) to describe the 
dynamics of solutions. It has been argued that Eq. ( [A7[ ) is inconsistent, in the sense 
that whereas the thermal fluctuations of the particles are taken into account, the thermal 
fluctuations of the solvent velocity field are not (the Oseen tensor gives the average solvent 
velocity as a function of the forces). In order to obtain the same level of description for both 



the particles and the solvent, Oono and Freed |TB[ introduced a set of coupled Langevin 
equations. However, it seems to us that the only difference between these Oono-Freed 
equations and Eq. ( |A^ ) lies in the neglect of the solvent inertia in the latter, for the 
following reason. It was shown in Ref. that if one derives the MSR functional from the 
Oono-Freed equations, the velocity and its conjugate appear quadratic in the Lagrangian, 
and so they can be integrated out explicitly. The resulting functional is identical to the 
one obtained from Eq. ([A^), provided that one takes the solvent density to be zero. Since 
the effects of solvent inertia on the dynamics of polymer solutions are negligible anyway, we 
conclude that Eq. (|A7|) is equivalent to the Oono-Freed equations. As an illustration of the 
irrelevancy of the solvent inertia, consider the Zimm model. The characteristic frequency lj 
of fluctuations with wave vector k is given by [^] = k-^Tk^ /Girr], where t] is the solvent 
viscosity. It follows from the Navier- Stokes equation that the solvent inertia is negligible if 
puj << rik"^, where p is the solvent density. Taking the viscosity and the density of water 
one finds that the effects of the solvent inertia are negligible on length scales L » lO^^^m 
. Nevertheless, if one is interested in the correlation and response functions of the solvent 
velocity field, the Oono-Freed equations are certainly usefuU [If. 
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FIGURE GAPTION 

The dynamic charge density correlation function g{q, s) as a function of the dimensionless 
time s := t/r for fixed dimensionless wave vector 0.1 < q := < 1, where ^ is the blob 
size, and r is the characteristic time scale at the blob level. 
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